Vibrational properties of single-wall nanotubes and monolayers of hexagonal BN 
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We report a detailed study of the vibrational properties of BN single-walled nanotubes and of 
the BN monolayer. Our results have been obtained from a well-established Tight-Binding model 
complemented with an electrostatic model to account for the long-range interactions arising from 
the polar nature of the material, and which are not included in the Tight-Binding model. Our study 
provides a wealth of data for the BN monolayer and nanotubes, such as phonon band structure, 
vibrational density of states, elastic constants, etc. For the nanotubes we obtain the behavior of the 
optically active modes as a function of the structural parameters, and we compare their frequencies 
with those derived from a zone-folding treatment applied to the phonon frequencies of the BN 
monolayer, finding general good agreement between the two. 

PACS numbers: 63.22.+m, 62.25.+g, 78.30.Na, 78.20.Bh 
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I. INTRODUCTION 

The discovery of the fullerenesS iruthe mid 1980's, and 
especially that of carbon nanotubeso in the early 1990's 
has fueled a frenzy of research activity in the field of 
carbon nanostructures, doubly motivated by their large 
potential for practical applications on one hand, and be- 
cause these systems constitute a new open field for basic 
research at the nanoscale on the other. Nanotubes can 
be single- walled or multi- walled, and are quasi- ID struc- 
tures having large aspect ratios (the quotient of length 
over diameter). Single- wall carbon nanotubes (SWNTs) 
can be either semiconducting or metallic, depending on 
their structural parameters, and can be chiral or achiral. 
They have large thermal conductivity ,|-aijd-have the high- 
est Young's modulus ever measurednoDu. Nanotubes 
can be filled with other elements, compounds or even 
fullerenes in their internal channel, forming structures 
similar to nanoscale coaxial cables. All these proper- 
ties confer nanotubes a large potential for technological 
applications, some of which have already been demon- 
strated. Nanotubes have been used to fabricate field 
emission devicesu, tips for scanning probe microscopy in- 
struments], constituents of nanoekctronic devicesH, gas 
sensorsEj, composite reinforcementEJ, lubrication, etc. It 
is therefore not surprising that these systems have been 
and continue to be the focus of .considerable interest: the 
large number of review articlestHI and monographslij de- 
voted to the different aspects of fullerenes and nanotubes 
in recent years provides further proof of this interest. 

One of the probes most frequently used to characterize 
carbon nanotube samples is vibrational spectroscopy, and 
specifically Raman spectroscopytLa. Raman spectroscopy 
offers experimentalists a rapid way of estimating the di- 
ameter distribution of tubes present in a sample, be- 
cause the low frequency radial breathing mode (RBM) of 



SWNTs has a frequency which is inversely proportional 
to the square of the nanotube diameter, independently 
of the nanotube structure, determined by the nanotube 
indices (n,m). Also, the high frequency G band modes re- 
sulting from atomic vibrations in the nanotube shell, can 
help to distinguish between metallic and semiconducting 
nanotubes in the sampletJ. Jorio and coworkerst£l have 
recently demonstrated the ability of Raman spectroscopy 
to provide full structural determinations of isolated car- 
bon SWNT's. The diverse aspects of phonons in carbon 
SWNT's and the use of Raman spectroscopy as a tool 
for their characterization-have been recently reviewed by 
Dresselhaus and Eklundta. 

Soon after the discovery of carbon nanotubes it was 
speculated that other materials could possibly form sim- 
ilar nanostructures, since there are several elements and 
many compounds which form layered structures bearing 
some resemblance to graphite. The most obvious can- 
didate was hexagonal BN (h-BN), which-, was predicted 
on the basis of theoretical calculationsEZltLa to be capable 
of forming nanotubes, a prediction which was later cor- 
roborated experimentally by the synthesis of such nan- 
otubesta. Today we know thatmapj/ other structures can 
form nanotubes (M0S2, WS2EI, BirB, . . . ); nevertheless it 
is still the case that carbon nanotubes continue to attract 
a larger interest, but certainly these other structures are 
interesting in their own right, and may be able to offer 
different possibilities for technological applications that 
carbon nanotabes cannot provide. Both multi-walled and 
single- walledE3 BN nanotubes can now be readily synthe- 
sized, and these tubes are uniformly insulating, tending 
to have a zig-zag structure. 

As other types of nanotubes become more common, 
many of the characterization tools and techniques ex- 
tensively used in the case of carbon nanotubes will un- 
doubtedly also find application in the study of these other 
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structures. In particular, vibrational spectroscopy, which 
has proved to be such a useful tool in the case of carbon 
nanotubes, is likely to prove useful also in these cases. 
In this context, it is interesting to consider the phonon 
properties of BN nanotubes from a theoretical point of 
view; their study will help to develop the characteriza- 
tion potential of spectroscopic techniques when they are 
applied to BN nanotubes. The aim of this paper is to 
provide such a theoretical study. We report extensive 
theoretical calculations of the vibrational properties of 
a h-BN layer and of a number of BN nanotubes having 
diameters in the range of 0.4 to 2 nm, including zig-zag, 
arm-chair and several chiral nanotubes. From our results 
for the flat sheet we have performed a zone-folding anal- 
ysis in order to predict the vibrational properties of the 
tubular structures from those of the flat layer, a tech- 
nique which has been frequently used in the case of car- 
bon nanotubes. The direct calculation of the vibrational 
properties of tubular structures allow us to compare the 
predictions of zone-folding with the actual results, and 
thus gauge the applicability of the zone-folding approach 
for BN nanotubes. Other authors-h.ave previously consid- 
ered the vibrationals and elasticE-3 pxaperties of isolated 
h-BN monolayers and of bulk h-BNc§E!± Some aspects 
of the elasticity of BN nanotubes have also been studied 
in references pi\ and |24|. However, the study of the vibra- 
tional properties in the case of BN nanotubes has not, to 
our knowledge, been addressed to date. 

The structure of the paper is as follows. In sec. [0| wc 
describe the model and calculation procedure used in our 
study; later in sec. Ill we discuss our findings, starting 



with a thorough description of the vibrational proper- 
ties of an isolated h-BN monolayer, followed by the re- 
sults obtained for the nanotubes, establishing a compar- 
ison between the results obtained from the zone-folding 
analysis and those obtained from direct calculation. Our 



summary and conclusions are discussed in sec. [V 



II. 



MODEL AND COMPUTATIONAL 
PROCEDURE 



All calculations have been perihrmed using the non- 
orthogonal Tiaht-Binding (TB)Ej parametrization of 
Widani et alB3. This is a parametrization following 
the scheme proposed by Porezag and coworkers^. This 
model employs a basis set consisting of one S and three P 
functions per atom, to represent valence states, and the 
resulting Hamiltonian and overlap matrix elements ex- 
tend up to a range of approximately 5 A. The core-core 
repulsion is modeled by means of a simple pair-repulsive 
potential. The original work of Porezag et alEi described 
a non self-consistent TB model; although an extension of 
the methodology has been proposed which incorporates 
some degree of self-coasistencyEJ, here we only have used 
the model of WidanycJ in its non self-consistent form. It 
is worth commenting that in spite of its ..simplicity this 
model has proved to be quite successfulEZrEa. 
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FIG. 1: Force constants derived from the TB model (symbols) 
and the typical behavior of the electrostatic interactions (ap- 
proximated in the plot by f(r) Zy/(er 3 )) (solid line), as a 
function of distance. The inset shows more clearly the behav- 
ior for large distances. 



The equilibrium structures of the tubes were obtained 
by careful minimization of the total energy with respect 
to both the atomic coordinates and the lattice constant 
along the tube axis. For the phonon calculations, wc 
then computed the force-constant, matrix in real space 
using a finite-difference approachlEl. We used atomic dis- 
placements of 0.02 A, and the force constants were taken 
as the average of the results obtained with positive and 
negative displacements, in order to eliminate anharmonic 
effects. Since all nitrogen and boron atoms are equiva- 
lent by symmetry in the tubes, we only calculated the 
force constants for one nitrogen and one boron atom in 
the supercell and generated the rest of the matrix using 
the symmetry operations. 

The force-constants have to be calculated between a 
given atom and all the other atoms in the system. How- 
ever, since the interactions in our TB model decrease 
rapidly with distance, only the elements with atoms suf- 
ficiently close have to be computed. To do so, we set up a 
supercell large enough that a sufficient number of neigh- 
bors is included. It must be kept in mind that, in the 
supercell geometry, a given atomic displacement in the 
central cell is always accompanied by the same displace- 
ment of all image atoms. We need, therefore, to use a su- 
percell such that the effect of the image displacements is 
negligible. For this purpose we have used supercells con- 
taining six unitcells in the case of (n,n) and (n,0) tubes, 
and two and one unitcells, respectively, in the case of the 
(10,5) and (10,7) tubes. The supercells contain hence a 
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minimum of 144 atoms for the (6,0) and (6,6) tubes, and 
a maximum of 384 atoms in the case of the (16,0) tube. 

The above described TB model does not by itself in- 
corporate any long-range electrostatic interactions. Such 
interactions are important when considering the vibra- 
tional properties of polar materials, because they will 
make a long-range contribution to the force constants. 
Therefore we have opted for correcting this shortcom- 
ing in the TB model in a physically sound but other- 
wise empirical way. The key observation in order to 
introduce the effects of the electrostatic interactions is 
that, when an atom i suffers a displacement u(z) from 
its equilibrium position, a net electric dipole of magni- 
tude p J = YV Zl^u 1 ^ appears associated with this atomic 
mo^mentlM Here Z % Vii is the Born effective charge ten- 
soiH of the atom i. As a consequence, if two atoms i and 
j are simultaneously displaced from their positions, be- 
sides the interaction energy given by our TB model, it is 
necessary to include the long-range interaction between 
the two electric dipoles generated. This gives rise-Jto a 
new term in the force-constants matrix of the formEj 

C™{vij) = (1) 

where the superindex LR in C~*(rjj) indicates that this 
is the long-range contribution to the force constants; 
f{rij) is a switching function (see below), e is the dielec- 
tric constant, is the vector going from atom i to atom 
j, and ry is the distance between these two atoms. The 
switching function is designed to come into play at dis- 
tances sufficiently large so as to not affect the TB model; 
it takes the form 

f(r) = l-e-^ 3 . (2) 

Following the usual approacbEi, it is also necessary to 
modify the on-site elements of the force matrix in order 
to satisfy the acoustic sum rule. Therefore, we take 

C^f(0n) = -E^(^)- ( 3 ) 

The Born effective charge tensor has been taken from 
recent ab initio density functional calculationscS Z± = 
-Zf=0.82 for out-of-plane, and Zf= -Z | j v =2.71 for- 
m-plane displacementsEJ. We have taken e = 4 and 
r c = 2.6 A, the distance of separation between second 
nearest neighbors in the h-BN layer. We have not at- 
tempted to fit these values in order to improve agreement 
of our results with experimental measurements. Rather, 
our approach has been to take reasonable values for them 
in order to estimate the effect of the polar interactions 
on the vibrational properties of the BN monolayer and 
in the nanotubes. Our chosen value for e is close to the 
experimental one for h-BN (4-5 depending on the direc- 
tionEaEa), and the value of r c used in Eq. (||) is moti- 
vated by the fact that force constants arising from the 



TB model alone decay very rapidly with distance, and 
they are already very small at this value of r c , as can be 
seen in Fig. [3. By choosing r c equal to the second nearest 
neighbor distance we ensure a smooth switching between 
the short-range TB model and the long-range electro- 
static one, with minimum interference between them. 

Once we have the force-constant matrix in real space, 
we can calculate the dynamical matrix in reciprocal 
space, and diagonalize it to obtain the phonon modes 
and frequencies. In the case of bulk polar systems the 
computation of the dynamical matrix in reciprocal space 
has to be made with some care. The reason is that, due 
to the long-range nature of the dipole-dipole interaction 
in Eq.(|^), a different limit for the dynamical matrix is 
obtained, as the phonon wavevector k approaches T, in 
the cases of longitudinal and transversal vibrations. In 
our case however, due to the reduced dimensionality of 
the systems considered (a monolayer and single-walled 
tubes), the Fourier transform of the dynamical matrix 
can be performed without further complications by sim- 
ply adding up the elements of the real-space matrix with 
the appropriate phase factors. Given its simple analyti- 
cal form, we can include a large number of neighbors (we 
tipically take all neighbors within a radius of ^200 A) in 
the summation involving the dipole-dipole interactions 
and guarantee a good convergence. The same limit is 
found for both polarizations when k — > T, and there are 
no splittings between in-plane longitudinal and transver- 
sal optical modes at T (see Appendix ). 

III. RESULTS AND DISCUSSION 

One of the aims of this work is to judge the validity 
of the zone-folding approach for obtaining phonon in- 
formation of tubular structures from the phonon band 
structure of the flat sheet, a procedure which has. been 
frequently used in the case of carbon nanotubesEJ. To 
this end we have performed a phonon analysis of the 
fiat sheet, and of several nanotubes (arm-chair, zig-zag 
and chiral) so as to compare the zone-folding results with 
those obtained directly from the tubular structures. In 
the following we will first present the results obtained for 
the flat sheet, followed by those for a selected set of nan- 
otubes. The comparison between these latter results and 
the zone-folding predictions will be then presented. 



A. Flat sheet 

Using the TB model discussed in section E3, with and 
without the dipolar interactions, with a supercell con- 
taining 192 atoms we have obtained the phonon disper- 
sion curves and vibrational density of states for a flat 
sheet of h-BN, plotted in Fig. |[ In order to check the 
convergence of the results with respect to the supercell 
size, we performed a test with a supercell containing 320 
atoms; the results obtained with the smaller cell are in- 
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FIG. 2: (a) Phonon band structure of the BN hexagonal monolayer (right) and the corresponding density of vibrational states 
(left); thick lines correspond to band structure and density of states calculated with electrostatic interactions included, while 
thin lines give results obtained only with the TB model, (b) Projected density of states on the B atoms; continuous line 
represents the density of states projected onto the plane of the monolayer, while the discontinuous line gives the density of 
states projected perpendicular to the plane of the monolayer, (c) As in (b) for the N atoms. 



distinguishable from those of the larger cell, and hence 
all subsequent analysis was performed with the 192 atom 
cell results. We also performed a more detailed analy- 
sis of the convergence of the phonon bands with respect 
to the lower threshold of the elements of the real-space 
force-constant matrix included in the construction of the 
dynamical matrix. We focus first on the convergence with 
respect to the TB force constants. We found that includ- 
ing all the couplings coming from the TB hamiltonian 
with absolute value larger than 1.2 x 10~ 2 eV/A was 
necessary to adequately converge all bands. This thresh- 
old corresponds to a spatial cutoff of 6.7 A, which in- 
cludes all atoms within the tenth nearest-neighbor shell. 
This large cutoff is required to correctly reproduce the 
quadratic behavior of the lowest acoustic band [labeled 
ZA in Fig.|| in the limit k — > T, with a zero sound veloc- 
ity (slope of the band) as k approaches I\ This quadratic 
behavior of the ZA band is a consequence of the fact that, 
at least to the lowest order in its amplitude, the strain 
energy created by this vibration is solely associated with 
the curvature that this out-of-plane bending mode in- 
duces in the layer. Optical bands typically converge much 
faster (ca. including interactions with neighbors within 
the third or fourth nearest-neighbor shell). It is worth 
mentioning in this respect that in the earlier theoretical 



work of Miyamoto et alE2 a supercell which only permit- 
ted the inclusion of second nearest neighbor interactions 
was used due to computational limitations (Miyamoto 
and coworkers used first-principles methods, which are 
significantly more costly, to evaluate the force-constant 
matrix). Although inclusion of up to second nearest 
neighbor interactions is sufficient to adequately converge 
most bands, it fails to reproduce the over-bending (see 
below) of the longitudinal-optical (LO) band, for which 
the maximum occurs not at T but approximately half 
way along the symmetry lines T — M and T — K, and is 
insufficient to reproduce the quadratic behavior of lowest 
acoustic band around T. 

Let us now comment briefly on the convergence with 
respect to the long-range electrostatic interactions. As 
already commented in the previous section, due to its 
simple analytical form we have used a large cut-off of 
200 A for the dipole-dipole interactions. However, we 
have checked that already with a cut-off of ^25 A the re- 
sults are well converged. Further reducing the interaction 
range starts to have an appreciable effect. Not only the 
position of the optical branches at and close to T start 
to shift from their converged values, but the finite slope 
at r that the LO mode must exhibit as a consequence of 
the long-range nature of the interactions also disappears 
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C m = C u 0.343 (0.322) 

Cxxyy = Cl2 0.119 

Cxyxy = C 66 0.107 (0.123) 
Y 0.302 
a 0.347 



TABLE I: Elastic constants and Young modulus Y, in units of 
TPa nm, and Poisson ratio a of the fiat BN hexagonal sheet, 
as calculated directly from the force-constant matrix given by 
our model. The results obtained from the sound velocities of 
the TA and LA bands (see Eq.(H)) are shown in brackets. 



(see Appendix ). 

As can be appreciated in Fig. due to the two- 
dimensional character of the system, there is a separa- 
tion between in-plane and out-of-plane modes, the latter 
being the lower frequency ones, with a maximum at T of 
719 cm -1 with dipolar interactions included (706 cm -1 
if these are not included). This frequency separation is 
to be expected, since in-plane modes excite both bond- 
stretchings (hard) and bond-bendings (soft), while out- 
of-plane modes result mostly in bond-bending motion 
and very little bond-stretching. The highest frequency 
bands, the LO and the transversal-optical (TO), have a 
frequency at T of 1446 cm -1 when the dipolar interac- 
tions are included (1479 cm -1 without). The maximum 
of the LO band occurs not at T, but at an intermediate 
point away from the Brillouin zone edges, as pointed out 
above. This is known as over-bending, and the presence 
of over-bending in the flat sheet phonon band structure 
can have interesting consequences for the phonon band 
structure of the nanotubes, as it may lead to the ap- 
pearance of new modes at T having higher frequencies 
in the tubular structures than in the flat sheetO. Note 
that including the dipolar interactions increases the over- 
bending, which is much smaller when the electrostatics 
are not appropriately accounted for. Another notewor- 
thy characteristic of the phonon band structure displayed 
in Fig. |2](a) is the splitting of approximately 250 cm -1 
existing between the ZO and ZA bands along M-K. In 
a homopolar hexagonal sheet such as graphene these two 
bands cross at KE3. 

From the slopes of the transversal-acoustic (TA) and 
longitudinal-acoustic (LA) bands in the limit k->Twc 
deduce a sound velocity of 13 Km/s and 21 Km/s, re- 
spectively. As already pointed out above, the ZA band 
exhibits a k 2 behavior in this limit, as expected, and 
therefore the sound velocity associated to this band is 
zero. Fitting the frequencies below 250 cm -1 to the ex- 
pression ui=2irv=5k 2 +~fk 3 , we can estimate a value of 
S w 1820 ± 60 cm- 1 ! 2 ^ 5.5±0.2 xlO^mV 1 . This 
parameter is interesting because it allows to estimate the 
energy necessary to roll up the BN sheet in order to form 
the nanotubes. Simple arguments show that the strain 
energy per atom can be approximated by E s t=C/r 2 with 
C=5 2 (mjv + ms)/4, being mpf and ijib the masses of 



the nitrogen and boron atoms, respectively, and r the 
tube radius. The value obtained for C in this way is 
1.91 eVA 2 , which is in reasonable agreement with the 
magnitude obtained from total energy calculations per- 
formed using both, a-iight-binding hamiltonian similar to 



the one utilized he 
tional calculational3'E 



and first-principles density func- 
A similar estimation in the case 
of grapheneEZI (6 ~ 6 xl0~ 7 m 2 s~ 1 ) leads to a value C« 
2.3 eVA 2 . Therefore, already looking at the phonon band 
structures it is possible to find an indication that the 
strain energy is smaller in the case of BN tubes than for 
carbon nanotubes, a result i}*daich is confirmed by more 
sophisticated calculationsMESEZIH 

Using the calculated force-constaate matrix we can di- 
rectly obtain the elastic constantsEa of the flat sheet, 
which are given in Table Q. In bulk three-dimensional sys- 
tems it is customary to quote elastic constants in units 
of pressure, due to the inverse equilibrium cell volume 
fartm- that appears in the definition of the clastic con- 
stantsHcj. However, in two-dimensional systems the 
definition of the cell volume is arbitrary, and it is there- 
fore more appropriate-,to use cell area in the definition of 
the elastic constantsEH; hence, we provide our results in 
units of pressure times length. Comparing with previous 
calculations in the literature, our result for C xxxx is in 
reasonable agreement with the value of 0.271 TPa nm ob- 
tained for the in-plane stiffnessEj by Kudin and cowork- 
ers using, density-functional theory and Gaussian type 
orbitalsE-i Our result is also in quite good agreement 
with the value of 0.309 TPa nm that can be deduced 
from the plane-wave density-functional calculation of the 
elastic constants of bulk h-BN reported by Ohba et al. in 
Ref. E6| (we use here the calculated interlayer distance of 
3.25 A to translate the data from the bulk to the mono- 
layer geometry). However, the comparison is poorer for 
the C xyxy elastic constant, for which the data of Ohba 
and coworkers translate into 0.055 TPa nm for the mono- 
layer (roughly half our value). The origin of this large 
discrepancy is unclear, although it might be related to 
the interactions with neighboring layers, which are not 
present in our case. 

For a hexagonal two-dimensional crystal only two elas- 
tic constants are requirecO, and due to the underlying 
crystal symmetry we obtain the following rule: 



= c 



xxyy 



-(- 2,C X y X y. 



(4) 



This provides an internal consistency check for our re- 
sults, which, as can be seen from the results in Table |, 
is satisfactorily obeyed [the numerical value obtained for 
Cxxxx using Cxxyy and C xyX y in Eq. (§ is within 3 % 
of its calculated value, given in Table Bj. A second con- 
sistency check is provided by the comparison of these 
values of the elastic constants with those obtained from 
the sound velocities reported above using the identities 



V LA 



Cxxxx 



„ 2D 
J TA 



Cxy 



(5) 
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where p is the surface mass density of the sheet. These 
data are also given in Table |. As can be seen there, the 
values obtained in these different fashions are in rather 
good agr-eement. From the elastic constants we can also 
calculatee3 the Poisson ratio, 



a 

Cx 



xxyy 



and Young's modulus, 



Y = 



^ xxxx ^ 'xxyy 
Cx.t.t.t. 



(6) 



(7) 



The values obtained for a and Y are also given in Table @. 
These values can be compared with those obtained pce.- 
viously for BN nanotubes using the same TB modclcH; 
the value of a obtained here is a little bit larger than 
that obtained for the nanotubes (0.263), but the agree- 
ment with the tube Young's modulus (0.298-0.310 TPa 
nm depending on the tube diameter) is nearly perfect. 
We also point out that the inclusion of the dipolar inter- 
actions does not affect the values of the sound velocities 
or elastic constants, as expected from the fact that they 
only influence the high frequency optical phonon bands, 
but not the acoustic ones, as can be seen in Fig. ^. 

At the zone edge the ZA band corresponds mostly to 
N atom out-of-plane displacements, while the ZO mode 
corresponds to B displacements. This is most clearly seen 
in Fig. ||(b) and (c), which shows the partial vibrational 
density of states (DOS) separated into its in-plane and 
out-of-plane atomic contributions. 

The frequency of the Raman and infrared-active (IR) 
modes mJaulk— h-BN has been investigated by several 
authorsHc3c3c3 obtaining similar results. The values 
reported . for . the high energy Raman-active E2 9 mode 
are 1366BN, 1367E3, and 1370El cm" 1 . For the IR 
modes the values given in the literature are: .1367E3 and 
1383±5El cm" 1 for the E lu (TQX mode, 16J0EI cm" 1 for 
the Ei„(LO) mode, and 770±3Band 783£2cm- 1 for the 
A2«(TO) mode. At this point it is important to note 
that, since our calculations refer to a free-standing (iso- 
lated) flat sheet of h-BN rather than to the bulk ma- 
terial, the comparison with these experimental observa- 
tions has to be made with some care. In analogy to 
the case of graphite, and assuming relatively weak inter- 
layer interactions, the frequencies of a single BN sheet 
should be closely related to those of the bulk. However, 
there is an important exception related to the role of the 
long-range electrostatic interactions. In a strictly two- 
dimensional (2D) system like the one treated here (i.e. a 
sheet of monatomic thickness), the electrostatic interac- 
tions do not result in a macroscopic term at T for in-plane 
modes, and hence we do not obtain any LO-TO splitting 
in that limit (see Appendix ). This striking difference 
is quite notorious, for example, comparing our phonon 
band structure in Fig. ^|(a) with the calculations for bulk 
h-BN diplayed in Fig. 3 of Ref. || The absence of LO- 
TO splitting for the systems of reduced dimensionality 



considered here turns out to be quite important when 
considering the applicability of the zone- folding approach 
to the vibrational properties of BN nanotubes, a point to 
which we will return below. 

Our results for the planar BN sheet have, therefore, 
to be compared with the TO frequencies of the bulk. 
We obtain for the highest modes at T a frequency of 
1446 cm -1 , to be compared against the empirical val- 
ues of the Fi2g and Ei„(TO) modes, and 719 cm -1 for 
the ZO mode, to compare with the measurements for the 
A2«(TO). The discrepancies are therefore smaller than 
10 %, which we regard as acceptable, given the sim- 
plicity of the TB model employed in this work. It is 
noteworthy that the effect of the inclusion of the long- 
range electrostatic interactions improves the comparison 
with experiment, correcting the under-estimation of the 
ZO(r) frequency and the over-estimation of the TO(T) 
and LO(r) value, as well as giving a more realistic over- 
bending of these modes along T — M and K F. We 
should emphasize here that no empirical information on 
vibrational or structural properties has been used-in the 
construction and parametrization of the modeled. As 
pointed out above, it is also important to take into ac- 
count that at least some of the discrepancy may not be 
attributed to the model, but to the geometry used in 
the supercell calculation. The presence of near-by sheets 
in the experiment is likely to soften somewhat the in- 
plane LO and TO modes, while hardening the ZO mode, 
which could account for some of the discrepancy between 
our results and the experimental ones. Discrepancies of 
similar size also occur between the first-principles values 
of Miyamoto et alt3 and the experimental values. The 
overall topology of our phonon band structure is in good 
agreement with the results of Miyamoto et a/., although 
we cannot make a direct comparison due to the unusual 
shape-pf the h-BN supercell used in that work. Kern 
et alEB have also performed first-principles density func- 
tional calculations of the phonon properties of both cubic 
and h-BN. Our phonon band structure is in reasonable 
agreement with that obtained by them, bearing in mind 
the different methodologies and geometries used in their 
work and ours. 

Rokuta et ai.@ have obtained the phonon spectra of a 
BN monolayer deposited over a series of metal surfaces 
(Ni, Pd and Pt) using high resolution electron energy loss 
spectroscopy. Again, there is overall good agreement be- 
tween the results reported by Rokuta and coworkers, and 
our own. Interestingly, they find no splitting at T be- 
tween the LO and TO modes when the BN monolayer is 
placed over a Ni(lll) surface, but some splitting is seen 
when the other metals are used. This result was con- 
sidered as quite surprising by these authors, especially 
taking into account that Ni(lll) was the sole substrate 
on which BN formed well ordered commensurate mono- 
layers, and it was suggested that this is because Ni can 
effectively screen the polarization field of the LO mode, 
while Pd and Pt cannot. However, we maintain that the 
polarization field only arises in 3D systems; in the 2D 
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FIG. 3: Comparison of the the phonon dispersion curves of 
the (10,10) tube (40 atoms in the unit cell and a diameter of 
13.9 A) calculated directly from the tube structure, and those 
obtained using the zone-folding approach from the results for 
the flat sheet shown in Fig. (0a). 



case the electrostatic interactions do not give rise to a 
macroscopic field for in-plane modes, i.e. the field has 
zero component at T, and hence there is no splitting of 
the TO and LO bands in that limit. However, the elec- 
trostatic interactions do give rise to non-zero components 
of the field in regions of the Brillouin zone away from T, 
and these result in a higher over-bending when the elec- 
trostatic interactions are included. Another consequence 
of the inclusion of the electrostatic interactions is the fact 
that the LO mode approaches T with a finite slope, which 
does not occur when the electrostatics are not included 
[see Fig. |(a)] . 

Prompted by this experimental study, we have consid- 
ered the effect of the presence of an ideal metal surface 
in the proximity of the BN layer on the phonon proper- 
ties of the latter. The metal surface is simply considered 
as a medium which generates dipolar images of the BN 
monolayer. The net result is to effectively shorten the 
range of the electrostatic interactions: the finite slope 
of the LO mode at T is lost, and although quadrupolar 
interactions persist in the BN-metal surface system, the 
dipolar interactions are screened out. 

Interestingly, Rokuta et alJiB report a degeneracy of 
the ZO and ZA modes at K. This is unexpected, because 
such modes should display a splitting related to the dif- 
ferent masses of B and N. The fact that no splitting is 
observed in the experiments is probably a result of a de- 
viation from planarity in the BN monolayer due to the 
presence of the metal surface. Another important point 
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FIG. 4: Same as Fig. | for the (14,0) tube (56 atoms in the 
unit cell and a diameter of 11.3 A). The vibrational density 
of states decomposed in the different directions of the atomic 
displacements is also shown. Solid lines for out-of-plane dis- 
placements and dashed lines for in-plane displacements. 



might be the quite different positions occupied by B and 
N atoms on the surface. In fact, Rokuta and coworkers 
reported a buckling of the BN monolayers on Ni(lll), 
where B atoms, which occupy fee sites, lie 0.2 A below 
the N atoms, which site directly on top of the Ni atoms 
in the surface. 



B. Nanotubes 

Let us now move on to discuss the results obtained 
for BN nanotubes. We have considered a series of (n,n) 
(arm-chair), (n,0) (zig-zag) and (n,m) (n ^ m =/= 0, chi- 
ral) nanotubes. For (n,0) nanotubes we have included 
all tubes with n from 6 to 16, i.e. diameters from 5 to 
14 A, approximately; while for (n,n) tubes we considered 
specifically the (6,6), (10,10) and (15,15) tubes, with di- 
ameters between 8 and 21 A. Only two chiral nanotubes 
were considered, since these structures usually have very 
large unit-cells, namely the (10,5) and (10,7) tubes, with 
diameters of 10.6 and 11.9 A, respectively. 

In Fig. | for the (10,10) tube, and Fig. §, for the (14,0), 
we compare the phonon dispersion curves calculated di- 
rectly from the tube structure, and those obtained for the 
same tubes employing the zone-folding approach from the 
phonon bands of the flat sheet discussed above. Firstly, 
let us remark that the zone-folding method cannot repro- 
duce certain low energy modes; in particular, it cannot 
give rise to a breathing mode, which in the flat sheet cor- 
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responds to a zero energy translation perpendicular to 
the plane. However, the high frequency modes are ex- 
pected (and indeed are) to be reasonably well described. 
In spite of the limitation regarding low frequency modes, 
the zone-folding approach has been extensively used for 
the interpretation of experimental results of carbon nan- 
otube vibrational modest3E3, partly due to its simplicity, 
and also to the fact that the phonon modes of graphite 
can be determined with very high accuracy. One of the 
aims of the present paper is to study if the zone-folding 
approximation can also be useful in the case of the BN 
nanotubes. 

Generally speaking, for experimentally relevant tube 
diameters, the zone-folding method reproduces quite well 
the overall phonon band structure. Deviations can be 
larger for the smaller nanotubes, where the effects of 
curvature are more noticeable. In particular, we shall 
see below how curvature effects influence the breathing 
mode. Another difference that can be appreciated be- 
tween the tube dispersion bands and those determined 
from zone-folding is a general softening of frequencies for 
the explicit calculation, where curvature effects are taken 
into account. The downward shift is especially clear at 
high frequencies, but is also noticeable at intermediate 
ones. This shift is not homogeneous, and it has also a 
non-trivial dependence on the wave vector and on the 
nature of the mode (radial or tangential). 

Also shown in Fig. ^ is the vibrational DOS for the 
(14,0) tube decomposed in the different directions of the 
atomic displacements. This curve resembles closely the 
result obtained for the BN monolayer. There is a clear 
energy separation between radial and tangential modes 
(parallel or perpendicular to the tube axis), the radial 
modes corresponding to frequencies below ~ 700 cm -1 . 
This corresponds to the already mentioned separation 
between out-of-plane and in-plane-modes in the case of 
the planar sheet, as could be expected, since the radial 
modes are roughly derived from the out-of-plane vibra- 
tions of the plane. 

There are four acoustic bands for the tubes. They cor- 
respond to a longitudinal mode (atomic displacements 
parallel to the tube axis), two degenerate transversal 
modes (atomic displacements perpendicular to the tube 
axis), and a torsion of the tube around its axis (the so- 
called twiston band which, in the limit of k=0, generates 
a rigid rotation of the tube around its axis). This last 
mode is characteristic of a one-dimensional system like 
the nanotubes, and no analog can be found for the bulk 
or the infinite planar sheet. According to Saito et aLcZI, 
and in contrast to the case of the flat sheet, all four bands 
would be required to approach T linearly. This relates to 
the fact that none of the acoustic bands of the nanotubes 
is solely derived from the quadratic out-of-plane acoustic 
band of the sheet: the transversal acoustic bands can be 
regarded as a combination of the TA in-plane and ZA 
out-of-plane modes of the planar sheet. However, it is 
interesting to point out here that the low energy behav- 
ior of the transveral acoustic bands in the nanotubes is 



still the subject of some controversy and some authors 
have recently proposed that these bands showJd exhibit 
a quadratic behaviour for small values of kE3't3. The sit- 
uation is more clear for the other acoustic bands since it 
is possible to establish a one-to-one correspondence be- 
tween the longitudinal and twiston bands of the tubes 
and the LA and TA bands of the planar sheet, respec- 
tively. In fact, the sound velocities of the latter bands 
are, within the precision of the calculations (~1 Km/s), 
independent of the tube radius and chirality, and almost 
identical (although somewhat smaller) to those of the cor- 
responding planar-sheet bands. Similar results, have been 
obtained for the case of the carbon nanotubesc3E3. This 
insensitivity to the tube structure and radius confirms 
the predictions of a simple continuum elasticity model 
of the tubesEil where the elastic constants are directly 
taken from the planar sheet, from which one obtains 
vjffi = w£5 an d v twiston = v ta- The independence of 
the mechanical properties on the chirality of the tubes 
is here a simple consequence of the isotropy of the un- 
derlying hexagonal structure of the BN planes. We have 
used this model to estimate the sound velocity of the 

transversal band, obtaining Vjf^ e — \J\ i>§^=9 Km/s. 
In the TB calculations, however, the transversal bands 
seem to be more sensitive to the numerical uncertain- 
ties than the other acoustic bands, being difficult to ex- 
tract an accurate slope. In fact, although the transversal 
bands apparently exhibit a linear dependence with the 
phonon wavevector (see figures ^ and |]), for very small 
values of k, and consequently for very small frequecen- 
cies, they seem to bend slightly . This would be roughly 
consistent with the proposedOEJ quadratic behavior, al- 
though problems of numerical accuracy cannot be ruled 
out. The clarification of this point will be the subject of 
future work. Here we have decided to make an estima- 
tion of the sound velocity of the transversal modes from 
the slope of these bands in the region where they still ex- 
hibit a clear linear behaviour. In this way we have found 
values ranging from 7 Km/s, for the narrower tubes, to 
8 Km/s, for those with larger diameters. This is again 
in reasonable agreement with the result that we obtained 
from continuum elasticity theory. 

In Fig. H we have plotted the frequencies of the op- 
tically active modes for the different tubes considered 
as a function of the tube radiusEa. The values plot- 
ted have been obtained including the long range elec- 
trostatic interactions, but the picture does not change 
significantly if they are not included^!. We have classi- 
fied the modes [according to the corresponding symme- 
try point groupE-3, which for zig-zag (n,0) nanotubes is 
Ginvi for arm-chair (n,n) nanotubes is Cmh-, and for chi- 
ral (n,m) nanotubes is Cn, where N=2(n 2 +m 2 +nm)/d^ 
with cW, being the greatest common divisor of 2n+m and 
2m+rla, i.e. N=70 for the (10,5) nanotube, and N=146 
for the (10,7) nanotube. The total number of Raman and 
IR active modes and their distribution over the frequency 
spectrum is very similar for both zig-zag and chiral nan- 
otubes. In the (n,0) nanotubes the rotationally invariant 
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FIG. 5: Active modes for the (n,0), chiral and (n,n) BN 
nanotubes considered in this work. Symbols have the same 
meaning for (n,0) and chiral tubes, except that for chiral tubes 
Al has to be read as A (see the text). See Ref. |5^ for the 
numerical values of the frequencies. 



modes can be classified according to how they transform 
under reflection through the symmetry plane a v , namely 
Ai and A2 modes; only the first are Raman and IR active, 
the A2 being inactive. The chiral tubes do not contain 
a cr v mirror plane, and therefore all six A modes are ac- 
tive. In total, the (n,0) nanotubes have eight modes with 
non-zero frequency which are both IR and Raman active 
(three belong to the Ai representation, and five to the 
Ei), and six modes (E2) which are only Raman active, 
while for the chiral tubes the count is nine IR and Ra- 
man modes Cfour A and five Ei), plus six (E 2 ) Raman 
active modesE-a. In the case of the (n,n) nanotubes, due 
to the existence of inversion symmetry in the point group 
the active modes can only be either IR or Raman active, 
having a total of four IR active (one A u and three Ei u ) 
and nine Raman (three A 9 , two Ei g and four E2 9 ) active 
modes. Comparing with carbon nanotubes, the bigger 
difference appears in the case of the (n,0) nanotubes, 
where the number of active modes is neatly double for 
BN nanotubes than for carbon nanotubesEi 

In Fig. ||(a) we have plotted the frequency of the 
breathing mode as a function of the tube radius. This 
mode is of importance in carbon nanotubes because it can 
be correlated with the tube radius and it is Raman active 
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FIG. 6: (a) Frequency of the radial breathing mode of BN 
nanotubes as a function of tube radius. The inset presents a 
log-log plot which clearly shows the r _1 behavior. Panel(b) 
shows the relative radial displacement of the B atoms and 
the N atoms in the radial breathing mode. The evolution 
of the buckling in the nanotube structure is plotted in the 
inset as a function of the radius. The meaning of the symbols 
is identical in all the panels and plots. See Ref. for the 
numerical values of the frequencies. 



for all tube chiralities. The same is true here. Although 
there are minor deviations at small radii (below 4 A) , the 
breathing mode frequency follows very closely a depen- 
dence on the inverse of the tube radius: vbm{t) = A/r. 
This is corroborated in the log- log plot shown in the inset 
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of Fig. ||(a), which fits very well to a linear behavior with 
a slope -0.97 ±0.01 (-1.09 ±0.03 including the tubes with 
r < 4 A). Leaving out the tubes with r < 4 A, a fit to 
this expression gives A = 1091 cm _1 A, with a standard 
deviation smaller than 3 cm -1 . This value is quite insen- 
sitive to t hc j. tu be structure, as occurs also with carbon 
nanotubesEZlE3. As can be seen in Fig. ^(a), the breath- 
ing mode exhibits a hardening for the smaller structures 
(r < 4 A) respect to what would be expected from the 
r _1 fit to the frequencies of the tubes with the larger 
diameters. If we include the frequencies of the smaller 
tubes, the quality of the fit decreases (mean deviation 
of 16 cm -1 ) and the value obtained for A increases to 
1160 cm _1 A, in accordance with the mentioned hard- 
ening. Interestingly, this behavior is the oppositte to 
what was observed for the carbon tubesEH, where there 
is a softening of the breathing mode for increasing curva- 
ture. Therefore, this effect must be related to an unique 
characteristic of the BN tubes as compare to the car- 
bon nanotubes: the buckled structure of its surface. For 
small and moderate radius BN nanotubes, the B atoms 
lie somewhat closer to the nanotube axis than the N 
atoms. This buckling is of the order of 0.1 A for the 
smallest rvatiotubes, but decreases rapidly with increas- 
ing radiuaHZl. Fig. ^(b) shows how the nanotube atoms 
are displaced in the breathing mode. The values plot- 
ted are the magnitude of the radial displacement of B 
atoms (5Rb) divided by that of the N atoms (6Rn)- For 
8Rb/5Rn = 1 the displacements are of equal magni- 
tude, and we would have a pure breathing mode, just as 
in a carbon nanotube. The observed behavior is slightly 
different: for tubes with r < 4 A the displacement of 
the N (outer) atoms is larger than that of the B (inner) 
atoms. This can be explained in a simple model as a 
consequence of the tendency to preserve in the pattern 
of displacements of this vibrational mode the bond angles 
of the buckled surface- 
As is well knownoc3C3, the constant A can be esti- 
mated from the elastic constants of the plane, 



(8) 



Using this equation we obtain the value A—11S5 cm _1 A, 
which is in reasonable agreement with the value found by 
the direct calculation. 

Finally, let us discuss the ability of the zone-folding 
approach to predict the frequencies of the active modes 
of BN nanotubes. In the figures |?] and || we plot the 
frequencies obtained from zone-folding side-by-side with 
those obtained directly from the nanotubes, in both cases 
including the electrostatic interactions. For clarity, this 
is done separately for intermediate frequencies (around 
700 cm" 1 ) and for high frequencies (above 1250 cm -1 ), 
and separately for (n, 0) and (n, n) nanotubes. Overall, 
we can say that zone-folding is capable of rather accu- 
rate estimations of the actual frequencies, particularly for 
nanotubes with radii larger than 4 A. In this case even 
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the ordering of the modes is accurately reproduced, both 
for (n, 0) and (n, n) tubes. The intermediate frequency 
modes are very well reproduced, but there is a slight over- 
estimation of the high-frequency modes, of the order of 
a few tens of wavenumbers. 

That the zone-folding approach is capable of providing 
accurate estimates of the BN nanotube phonon frequen- 
cies rests on the absence of a k = (macroscopic) dipolar 
field in the BN monolayer. The reduced dimensionality 
of the system prevents the build up of a k = com- 
ponent of the electrostatic field, which thus only has a 
finite range. Therefore, no splitting of modes is observed 
at T; the only effect of the electrostatic interactions is to 
shift slightly the high frequency modes and to confer a 
finite slope to the LO mode in the k — > T limit, modifi- 
cations which do not prevent zone-folding from providing 
accurate estimates of the phonon modes of isolated nan- 
otubes. A subtle point here is that the applicability of 
zone-folding in the case of BN also depends on the de- 
gree of transferability of the dielectric constant e and the 
Born effective charges [see Eq.(|])] between the monolayer 
and nanotube geometries. For small nanotube radii there 
may be noticeable changes in the dielectric constant, but 
this is unlikely to affect the application of zone-folding to 
all but the narrowest nanotubes. 



IV. SUMMARY AND CONCLUSIONS 

In this work we have presented a detailed analysis of 
the vibrational properties of a BN monolayer, and of 
single- walled BN nanotubes of different diameters, in- 
cluding (n,0), (n,n) and two chiral nanotubes. Our cal- 
culations are based on a non-orthogonal TB model com- 
plemented with long-range electrostatic interactions, not 
included in the TB model. The incorporation of these 
electrostatic interactions corrects to some extent the de- 
ficiencies of the TB model at intermediate and high fre- 
quencies in the BN layer. Using a simple model for a 
metal surface, we have analyzed the influence of placing 
a BN monolayer over such a surface on the phonon bands 
of the monolayer. The effect appears to be small and con- 
fined only to the high frequency range, above 1250 cm . 
Generally, our calculated phonon frequencies, sound ve- 
locities and elastic constants are in good agreement with 
previous experimental and theoretical results. Concern- 
ing the nanotubes, we have obtained phonon band struc- 
tures and have analyzed the optically active modes (IR 
and Raman) as a function of the nanotube structure and 
radius. We have also performed a detailed comparison 
between the predictions of the zone-folding approach and 
the actual frequencies of optically active modes obtained 
directly from the nanotubes, finding rather good agree- 
ment between the two for all except the narrower tubes, 
having radii below 4 A. High frequency mode frequen- 
cies are systematically overestimated by about a few tens 
of wavenumbers, but intermediate mode frequencies are 
given with high accuracy. Therefore, zone-folding should 



also be useful in the context of BN nanotubes, as it has 
been in spectroscopic studies of carbon nanotubes. We 
also report on the behavior of the frequency and the pat- 
tern of atomic displacements of the radial breathing mode 
as a function of the tube radius. 
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APPENDIX: ABSENCE OF LO-TO SPLITTING 
FOR A BN MONOLAYER 

In this appendix we briefly discuss the role of long- 
range electrostatic interactions in a monolayer, compar- 
ing the results with those found for bulk systems. Let 
us consider a thin slab of polar material parallel to the 
xy plane. We focus first in the in-plane optical modes, 
for which the atomic displacements u are parallel to the 
xy plane. Associated with this phonon there is an elec- 
tric polarization field P = j^(Zu)f(z) cos(k • r) and a 
charge density pp = — V • P, where Z is the Born effec- 
tive charge tensor within the unit cell, A is the area of 
the unit cell, and f(z) is a function related to the profile 
(number of layers) of the slab which can be approximated 
by S(z) for a monolayer. We stress that k, u, and P are 
all parallel to the slab for pure in-plane modes. For the 
particular case of a system with hexagonal symmetry, like 
our BN sheets, the in-plane properties are isotropic, and 
Z can be regarded as a scalar. Although this is not crucial 
for our arguments, we will make this simplification in the 
following. The charge density p p is proportional to the 
product k • u , and no macroscopic charge accumulations 
appear due to the TO modes. However, a macroscopic 
electric field E mac can appear associated with the LO 
modes. The correct limit when k — > is then calculated 
by solving Poisson's equation for p p . The components of 
the electric field in the directions parallel to the slab are 

E^ac = k*B(z, k) cos(k • r), (A.l) 

where k is the modulus of k, and B(z, k) is given by 

f°° f(ak)e~ i( - akz t 
B(z,k)= da m ^ 2 , (A.2) 

f(q) being the Fourier transform of f(z). Once the elec- 
tric field has been calculated, the corresponding mod- 
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ification of the force-constants matrix can be obtained 
taken into account that the force acting on an atom is 
given by F = Z E OTO( £3 . In the case of the monolayer 
we take j[z) — S(z), then B(z = 0, k) = 1/2 and the 
correction of the force-contants matrix for the in-plane 
modes becomes 



dj — 



2tt Z 2 k l W 



A 



(A.3) 



This correction goes to zero when k — > irrespective of 
the polarization of the mode and, therefore, both LO and 
TO in-plane modes have the same frequency. This is in 
agreement with the results ob tain ed for slabs of ionic ma- 
terials using different models 5 Tl 58 !: pure in-plane optical 
modes have frequencies close to those of the TO phonons 
in the bulk. It can also be noticed the linear behavior 
of the LO mode as it approaches V. We move now to 
analyze the bulk limit. By taking f{z) = l/L, being 
L the distance between atomic layers, we can study the 
limit for bulk phonons with k z = (k z is the wavevector 
perpendicular to the atomic layers in the bulk). B(z, k) 
becomes now l/(Lfc), and the term in the force^constants 
matrix takes the usual form for bulk systemsE3 
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which has a different value for LO and TO modes. For the 
intermediate case of a slab with finite thickness we find 
bulk behavior for phonons with sufficiently large k, while 
the behavior of the monolayer is recover when k goes to 
zero. The slope of the LO modes at F is proportional to 
the number of layers in the slab. 



For modes with polarization vectors having compo- 
nents along the normal to the slab (out-of-plane com- 
ponents within the nomenclature used throughout this 
paper) things can be quite different. For these modes, P 
has a component along z and, when k = 0, the charge 
distribution p p associated to them can be rationalized as 
a number of infinite planes with charges of alternating 
signs along z. There are also surface charges given by 
P ■ n = P z . All these charged planes give rise to elec- 
tric fields which do not decay along the entire thickness 
of the slab. In fact, those modes with atomic displace- 
ments perpendicular to the slab, although are transversal 
by definition (since the only phonon wavevector that can 
be safely defined is parallel to the slab), are known to 
have frequencies similar to thosq— Qf-,thc corresponding 
longitudinal phonons in the bulkfTEj. For this reason, 
Fuchs and KliewertiO and LucasE3 arrived to the conclu- 
sion that the phonon-frequency distribution of an ionic 
crystal slab converges very rapidly to that of bulk. How- 
ever, for anisotropic systems like h-BN, where the fre- 
quency range of the in-plane and out-of-plane vibrations 
is very different, the conclusion might be different. 



In the case of a free monolayer the separation between 
in-plane and out-of-plane modes is exact for an arbitrary 
wavevector k. There is only one out-of-plane optical 
mode (ZO) which, although its frequency is slightly mod- 
ified by the long-range interactions, cannot exhibit any 
splitting since it is a non-degenerate mode. The LO and 
TO in-plane modes do not exhibit splitting as explained 
above. 
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